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ABSTRACT 



Context. Cosmological A'-body simulations have revealed many empirical relationships of dark matter halos, yet the physical origin 
of these halo properties still remains unclear On the other hand, the attempts to establish the statistical mechanics for self-gravitating 
systems have encountered many formal difficulties, and little progress has been made for about fifty years. 

Aims. The aim of this work is to strengthen the validity of the statistical-mechanical approach we have proposed previously to explain 
the dark matter halo properties. 

Methods. By introducing an effective pressure instead of the radial pressure to construct the specific entropy, we use the entropy 
principle and proceed in a similar way as previously to obtain an entropy stationary equation. 

Results. An equation of state for equilibrated dark halos is derived from this entropy stationary equation, by which the dark halo 
density profiles with finite mass can be obtained. We also derive the anisotropy parameter and pseudo-phase-space density profile. All 
these predictions agree well with numerical simulations in the outer regions of dark halos. 

Conclusions. Our work provides further support to the idea that statistical mechanics for self-gravitating systems is a viable tool for 
investigation. 
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1. Introduction 

The investigation of structures and properties of dark matter ha- 
los is one of the important open questions in modern cosmology. 
Cosmological A^-body simulations have revealed many empiri- 
cal r elationships concerning dark halos, such as the densi ty pro- 
file (lNavarroetalJll997l: iMoore et alj[l999l: lEinastd[T965.) . the 
velocity dispersion, the anis otropy parameter, and the pseudo- 
phas e-space density profile dNavarro et alj|2010l; iLudlow et alj 
12010) . Yet the physical origin of these halo properties still re- 
mains unclear 

For about fifty years, many authors have attempted 
statistical-mechanical approaches to investigate the properties 
of self-gravitating syste r ns (to name a few, Antonov 196^ 
Lvnden-Belll 119671: [Shul [19781: iTremaine et ah 1986: Sridhai) 



1987t IStiavelli & BertinI 19871: IWhite & Naravanlll987l: iTsalh: 



1988; 'Sperael & Hernauist"l992; 'Soker"l996'; Kull et al "1997'; 
^ Nakamura 2000; Chavanis 2002; Hi orth & Williams 2010). 
$-H However, these attempts have encounter ed many f ormal diffi- 
5^ , culties, and little prog ress has been made (iMo et al-lfeOlQ) . 

In a recent work dHe & Kang||2010l ), we employed a phe- 
nomenological en t ropy form of ideal gas, first proposed by 
IWhite & NaravanI (Il987h . to revisit this question. First, sub- 
ject to the usual mass- and energy-conservation constraints, 
we calculated the first-order variation of the entropy, and ob- 
tained an entropy stationary equation. Then, incorporated in the 
Jeans equation, and by specifying some functional form for the 
anisotropy parameter, we solved the two equations numerically 
and demonstrated that the anisotropy parameter plays an impor- 
tant role in attaining a density profile that is finite in mass, en- 
ergy, and spatial extent. If incorporated again with some empir- 
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ical density profile from simulations, the theoretical predictions 
of the anisotropy parameter and the radial pseudo-phase-space 
density in the outer regions of the dark halos agree very well 
with the simulation data. Finally, we calculated the second-order 
variation, which reveals the seemingly paradoxical but actually 
complementary consequence that the equilibrium state of self- 
gravitating systems is the global minimum entropy state for the 
whole system, but simultaneously the local maximum entropy 
state for every and any finite volume element of the system. Our 
findings indicate that the statistical-mechanical approach should 
be a viable tool to account for these empirical relationships of 
dark halos, and may provide crucial clues to the development of 
the statistical mechanics of self-gravitating systems as weU as 
other long-range interaction systems. 

Despite its great success, we should point out that the spe- 
cific entropy in that work is not defined in a consistent way. 
On the one hand, as addressed above, the anisotropy parame- 
ter is of great significance to get a finite density profile, whereas 
on the other hand, we had to use only the radial pressure (or 
velocity dispersion) for the entropy definition. Any combina- 
tions of other directions' velocity dispersions that included the 
anisotropy parameter, failed to produce correct equations, and 
hence were rejected. Besides, another shortcoming of that work 
should be pointed out, which is the agreements between all pre- 
dictions and the simulations are restricted only to a narrow radial 
range of r/r_2 > 2. 

The aim of the current work is to remedy these defects. With 
the analogy of the isotropic Jeans equation, we introduce the ef- 
fective pressure instead of the radial pressure to define the en- 
tropy. In a similar procedure, we derive a new entropy stationary 
equation, which can be easily integrated to attain the equation 
of state of the equilibrated dark halos, which iHe & KangI (1201 Oh 
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failed to obtain with the radial pressure. It is this state equation 
that we employ to derive the halo density profile. Additionally, 
with this improved treatment, the radial range of the agreements 
between the predictions and simulations is increased by one or- 
der of magnitude. 

This paper is organized as follows. In Sect. |2]we define the 
effective pressure and derive the equation of state for the self- 
gravitating dark matter halos. In Sect. [3] we solve this equation 
approximately, and compare all results with the simulations of 
dark halos. We present the summary and conclusion in Sect.|4] 

2. Basic theory and formulae 

2.1. Isotropic Jeans equation and effective pressure 

After experienced the violent relaxation process (iLvnden-Belll 
Il967h . the self-gravitating system settles into a virial equilib- 
rium state, which is usually called the quasi-stationary state 
dCampa et al.ll2009l) . which is described by the Jeans equation 
in spherical coordinate^] as dBinnev & Tremain3l2008h : 



d(pcrj) per? 

+ 2/3— 
dr r 



dO 
'd7 



GM 



(1) 



where M(r) - J^^ Anp-pdr, O is the gravitational potential, yS - 

1 - crj jlcr^ is the anisotropy parameter, and cr? and cr? are the 
tangential and radial velocity dispersion, respectively. 

ForyS = 0, i.e. the isotropic case, the Jeans equation reduces 
to the following. 



dr 



d<D 

'd7' 



(2) 



where - pcr^ is the radial pressure. Equation (|2| resem- 
bles the hydrostatic equilibrium equation of an ordinary fluid, 
Vp = -pV<^, in that the radial pressure gradient, the buoyancy, 
resists the gravitational force of the whole system. Enlightened 
by this similarity, we generalize Eq. (|2]i to the anisotropic case 
by defining the effective pressure P through: 



dr dr r 



Then, the Jeans equation, Eq. ([TJ, can be re-written in terms of 
the effective pressure P as 



dP 
d7 



-p 



dO 
d7 



GM 



(4) 



Below we will use the eflFective pressure to construct the specific 
entropy form. 

2.2. Equation of state 

As mentioned in the introduction, in lHe&Kang|(l2OT0l) we em- 
ployed the entropy principle and derived an entropy stationary 
equation by introducing the specific entropy. 



This entropy form was first used by IWhite & NaravanI (Il987h 
for the isotropic self-bounded gas sphere, but we assumed it can 
also be applied to the velocity-anisotropic case. This is not a self- 
consistent treatment, and also it's hard to generalize to include 
the anisotropy parameter. 



' Thirougiiout tiiis work, splierical symmetry is always assumed for 
tire equilibrated dark matter halos. 



Based on the above similarity between Eqs. (|2]i and (|4]l, we 
replace pr in Eq. dU by the effective pressure P to make the 
entropy more appropriate for the anisotropic case, so that the 
total entropy is 



5t = r 47Tr^ps'dr = f 
Jo Jo 



4nr^p\n{P^'^p-^'^)dr, 



(6) 



where s' = ln(P^^^p"^^^). Subject to the constraints of mass and 
energy conservation, we calculate the firs t-order varia t ion of the 
entropy in a similar process as the one in iHe & Kand (l2010h . to 
obtain the new entropy stationary equation as 



3dlnP 5dlnp_ ^dP 
2 dr 2 dr pdr' 



(7) 



where A is the Lagrangian multiplier corresponding to the energy 
conservation. The left-hand side o f this equation is formally the 
same as that of iHe & Kand ( |2010|) . except that p,- is replaced by 
the effective pressure P. 

Equation O can be readily transformed into 



— = -- 4- -/I 

dP 5P 5 ' 

and can be directly solved as 

p = AP + fiP^'\ 



(8) 



(9) 



where is an integration constant, /i and A are related to the total 
mass and energy of the dark halo, that is, they can be specified 
by the total mass M, and total energy E, as fi - fiiM, E) and 
A - A{M, E), and vice versa. 

Equation ^ describes the relationship between p and P, 
which is exactly the equation of state of the equilibrated sys- 
tem, but P is also related to p through the differential equation, 
Eq. dU, thus Eq. (|9]l cannot directly provide us with the dark halo 
density profile. Below we will explore how to obtain the density 
profile from this state equation. 



(3) 3. Results 



From Eq. (|9]l we can see that p scales with P as p ~ AP at the 
center of the dark halo, but p ~ jjP^^^ at the outskirts, since 
both p and P are large at the center, but small at the outskirts of 
the dark halos. This suggests that, it would be advantageous to 
analyze the approximate solutions of these two cases, before we 
discuss the general solution of Eq. (|9]i- 

3.1. Ai^O,fi^O 

As explained above, this case is the approximation of Eq. (|9]l 
at the halo center If A is renamed as A = 1 /cr^ (we will see A 
must be positive in the following^ then with p - 0, Eq. ^ is 
re-expressed as 



(5) P^pcri 



(10) 



Incorporating Eq. (HJl, and differentiating both sides of the equa- 
tion with respect to r, we have 



a" 



,dp _ dO 
^d7 " ~^ d7' 



(11) 



- cr^ has the same dimension as the velocity dispersions cr^, or <r^, 
but they are completely different in physical meanings. 



Kang & He: explanation of dark matter halo properties 



3 




Fig. 1. Isothermal density profiles. The inner slope of the den- 
sity profile is dependent on different A. The Einasto profile, 
ln(p/p -2) = -2/a'((r/r_2)"-l), with a k 0.17 (see iNavarro et alj 
I20ld) . is also indicated for comparison. 



with the solution as <1) = -cr^ In p, w hich resembles the equation 
that the isothermal gas satisfies (see iBinnev & Tremaiii a 12001 
p. 303). Differentiating Eq. (fTTI) again, we have 



y2 



dlnp 
~d7 



) = -AnGAp. 



(12) 



We show the solutions in Fig. [T] with the case of A - 2 corre- 
sponding to the singular isothermal spher^E p ~ r^^. We can 
also see that A with A <2 bends the singular isothermal solution 
toward a centrally-cored density profile, and the size of the core 
depends on the value of A, but the slope at the outskirts of the 
dark halos remains unchanged. 

3.2. A = 0,1^^0 

In this case, the equation of state, Eq. (|9]l, is reduced to 



P = ti-'''p"\ 



(13) 



which resembles a polytropic process of a common gas. 
Differentiating this equation with respect to r, and incorporat- 
ing Eq. dUi, we get the following differential equation 



dr \ dr^ 



IQnG 



(14) 



As analyzed previously, this equation is just an approximation 
of Eq. (|9]l at the outskirts of the dark halos. As r reaches zero, 
we can see that the right-hand side of the equation also reaches 
zero, so that we can approximate the solution with a power-law 
form, p ~ r", with the power-index n - Q, -1.5. Then starting 
from the two power-law inner solutions at sufficient small r, we 
numerically evaluate Eq. ( fT4l l with different /i, and show all re- 
sults in Fig.|2] Since Eq. (fT4l l is only valid at the outskirts of dark 
halos, we can see that regardless of the non-uniqueness of the in- 
ner slopes at small radii, all resulting density profiles behave as 
if they were truncated by the yU-term of this equation, so that both 
the mass and the energy of the dark halos are not infinite. 



^ For the characteristic density p_2 and scale r_2, we set 4nG = p_2 
r_2 = 1. 




r/r_2 

Fig. 2. Polytropic density profiles. Panels (a) and (b) correspond 
to the two approximate inner power-law solutions, p oc r", with 
the power index as n - - 1 .5 and n -0, respectively. All density 
profiles are self-truncated, with the truncation radii varying with 
different jj. 
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Fig. 3. (a) The density profiles and (b) the logarithmic slopes 
corresponding to the density profiles in Panel (a). The profile 
and the slope derived from Eq. (flST l are indicated as "solution" 
in the figure, with A - 2.21 , = 0.36. The major fitting formulae 
of the density profile, NFW, Moore, and Einasto profile, are also 
shown for comparison. 
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3.3. Ai^O,j^i^O 

With the heuristic results from the above two subsections, we 
explore the complete solution of this general case. It may not 
be easy to solve Eqs. (|4]i and (|9]l exactly, so we try to find the 
approximate solution to the density profile. We treat Eq. ( fTOl l as 
a first approximation and substitute it into the derivative of the 
state equation, Eq. (|9]l, then we have 



dp dO 1 



3fi 



(15) 



1/cr^. We numer- 



in which again cr^ is the renaming of /I, as /I - ^ , j. 
ically evaluate this equation, and the resulting density profile is 
exhibited in Fig. [3^. We also indicate in Fig.[3]3 the logarithmic 
slope of the density profile, i.e., y = -dlnp/dlnr. In the inner 
region p indeed behaves as the isothermal solution of Eq. ( fTSl ). 
In the outer region, the density profile does not decline as fast as 
the polytropic solution of Eq. ( fT4l l because we only consider the 
first approximation, but it is still obvious that the whole density 
profile is roughly the superposition of the above two subsections' 
solutions. 

We can see that with an appropriate value of A and fj., both 
the resulting density profile and its log-slope agree very well 
with the fitting formulae at the radii r/r_2 > 0.2, the right part of 
the vertical dotted line in Fig. [3] At small radii with r/r_2 < 0.2, 
the solution is qualitatively accepted, but the agreement is not as 
good as the case at large radii. 

3.4. crl,l3, andp/crj profiles 

If the density profile p(r) is determined, we can obtain the ef- 
fective pressure P by integrating Eq. (|4]i. Numerical simulations 
indi cate that the anisotro py parameter, is usually a small num- 
ber (Na varro et al.ll20lb ). From Eq. Q we can see that when yS 
is small, we can approximate the effective pressure P as 



by which cr^. can be expressed as cr? 
From Eqs. ([T]l and ( fTSl l we obtain 



Pip. 



dp 



dr 



(1 -H3//o-^^^/5p2/5) dr' 



(16) 



(17) 



Then if cr?, p, and hence the logarithmic slope of the density 
profile y = -d In p/d In r are given, (3 can be expressed as 



2y6 = 



.6/5 \ 



2/5 



diner? 

y - -—^ ^ a(r)y - b(r), (18) 



din r 



where fl(r) = 1 - cr2/cr?(l -h 3yUcr775p2/5), b(r) = dlncr?/dlnr. 

We present a simple analysis of what Eq. (fTsT l implies. In the 
inner region of dark halos p is large, the ter m including n can be 
neglected, and from the simulation results (iNavarro et al.ll2010h 
we know both a(r) and b(r) vary slowly with r, so it seems that 
there is an approximately but not strictly linear relationship be- 
tween /? and y. The existence of this li near relation has been ob- 
served in both Hansen & Moord (1200 6) and the inner region of 
iNavarro et al. l (l2oToir Because this analysis will not be valid for 
regions at larger radii, the linear relationship should not e xist for 
the wh ole r egion. Indeed, the lat est simulations by Navai TO et al.l 
( 120101) and lLudlow etal] ( 1201 Ol) show that /3(r) is not the mono- 
tonic function of r, that is, haloes are nearly isotropic near the 




Fig. 4. The radial velocity dispersion cr?(r), the anisotropy 
parameter /3{r), and the pseudo-phase-space density Q{r) are 
shown in panels (a), (b), (c), respectively. The results based on 
the resulting density profile in Sect. 13.31 are indicated as "solu- 
tion" in the figure. The results based on the Einasto profile are 
also shown for comparison. The vertical scales of panels (a) and 
(c) are arbitrary. 



center, radially biased to the maximum at some radius and ap- 
proximately isotropic agai n in the outskirts, which i s against the 
earlier monotonic result of'Hans en & Moord (120061) . 

In the inner region, since a < 1 and b > 0, so from Eq. ( fTSl ) 
we must have y > 2/3, and in the outer region, since y > 2 and 
y6 < 1, we always have y > 2/3. So we have -v > 2/8 in the whole 
regi on, which is consis t ent wi th the result of lAn & Evans! (l2005h 
and lCiotti & Morgantil ( l2010 l). 

Next, we analyze the behavior of the pseudo-phase-space 



density Q(r) = p/crj (or p/cr^^^). From Eq. ( fTSl l we have 



3 diner? 
2 dlnr 



3b 
Y 



2/3 + b 



3b 

-■ (19) 



ding _ dlnp 
d In r d In r 

In the inner region we know from'Navarro et al.' (2010^ that /3, a, 
and b all vary slowly with r, so that d In Qjd In r roughly seems 
to be a constant. If we take lHansen & Moord (!2006l) 's result as 
a reference, a » 0.4, b ^ 0.32, and assume /3 ^ 0.12 on the 
average, then dlnQ/dlnr ^ -1.88, very close to the s pherical 
secondary-infall similarity solution ( Bertsc hingeil 1 985h . 

With the density profile derived in the previous subsection 
(see Fig. [3^), and by Eq. ( fT6b . we calculate the approximate ra- 
dial velocity dispersion, cr?, which is shown in Fig.|4^. The ve- 
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locity dispersion calculated with Einasto profile is also indicated 
in the figure for comparison. We can see that both the two re- 
sults roughly agree with the simulation results of Navarr oet alj 
(12010) when r/r_2 > 0.2, as is consistent with the accuracy of 
the resulting density profile of Sect. 13.31 The discrepancy is sig- 
nificant when rj r_2 < 0.2, which may suggest that there may be 
still some crucial physics that we have failed to capture in the 
current formulation. 

Given the density profile p{r), and the velocity dispersion 
o-^(r), and with the two parameters A (or cr^) and fi being the 
same as we computed the density profile, we can also obtain 
j3{r) via Eq. ( fTSl l as well as the pseudo-phase-space density 
Q{r). The results are exhibited in Figs.|4]3 and|4];, respectively. 
The corresponding results derived by using the Einasto pro- 
file are also indicated for comparison. From Fig. |4j3, we notice 
that the non-monotonicity of /3{r) is once again reproduced in 
ou r current work , which agr ees well with t he si mulation results 
of iNavarro etal] (1201 Q) and lLudlow etal] JIoTo,) . However, we 
find that both the position and the height of the peak of the 
jS(r) profile deviate slightly from the simulations. We speculate 
that this discrepancy may be originated in the approximation we 
adopted for Eqs. ( fTSb - (T% . The result based on the Einasto 
profile cannot be accepted for r/r_2 < 0.2, as is consistent with 
the previous results for both p and crj. 

As expected, the resulting Q{r) profile closely follows a 
power law, r near the center of the dark halo. However, 
the predictions do show a manifest curve-up deviation from 
this power law at the outskirts of d ark halos, whi c h has been 
observed by the simulation result of iLudlow et al.l ( 20101). and 



has al ready been reproduced by our previous study (iHe & Kand 
120101) . 

4. Summary and conclusion 

Cosmological A^-body simulations have revealed many empir- 
ical relationships concerning dark halos. Up to now, however, 
there are no robust physical origins for these empirical rela- 
tionships. In this work, we aim to provide a unified statistical- 
mechanical explanation to the relationships concerning the mat- 
ter density, anisotropy parameter and the pseudo-phase-space 
de nsity profile . The m ain steps of this work are similar to those 
of IHe & KangI (l20Tol) . i.e., by using the entropy principle with 
mass- and energy-conservation as constraints to de rive the en- 
tropy stationary equation, Eq. (|7]i. Different from iHe & KangI 
however, based on the analogy with the isotropic Jeans 
equation, we phenomenologically introduce the efifective pres- 
sure P instead of the radial pressure p,- to construct the specific 
entropy. The effective pressure is a key quantity in this work and 
proves to be more powerful than pr- 

The entropy stationary equation can be easily integrated to 
attain Eq. (|9]l, the equation of state of the equilibrated dark halos, 
which He & Kang (2010) failed to attain with the radial pressure 
Pi-. It is this equation of state that we employ to derive the halo 
density profile p(r). Notice that this density profile has nothing 
to do with /3. 

First, we study the approximate solutions of Eq. (|9]l in two 
asymptotic cases, i.e. the inner and outer region of the dark ha- 
los. The inner solution is the isothermal density profile, with the 
inner logarithmic slope varying with different parameter A. The 
outer solution is the polytropic density profile, self-truncated at 
sufficiently large radii with a different parameter /v, so that the 
density profile has finite mass, energy and extent. The approxi- 
mate complete solution can be roughly regarded as the superpo- 
sition of the two asymptotic solutions. By choosing appropriate 



values of A and /i, we can reproduce both the density profile and 
its logarithmic slope quite well in the region r/r_2 > 0.2. In the 
inner region of r/ r_2 < 0.2, the solution seems to be qualitatively 
accepted, but the agreement between the solution and the fitting 
formulae is not quantitatively well. We speculate that there may 
be still some crucial physics that we have failed to capture in the 
current formulation. 

Then, with the resulting density profile, we compute the ra- 
dial velocity dispersion cr^(r), and then /3(r) and the pseudo- 
phase-space density Q(r) profile. We find that all these predic- 
tions agree quite well with the corresp ond ing results of the lat- 
est simulations of Navarro et al. I (120 10) and lLudlow etall (1201 Oh 
in the region of r/r_2 > 0.2. Despite the slight discrepancies 
from the simulation results in the inner regions, these agreements 
at the ou ter regions of dark halos strengthen the conclusions 
of iHe & Kang (2010), indicating once again the great success 
of the statistical-mechanical approaches for the self-gravitating 
systems. 
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